NAVW.  P0S1GRA0UATE  SCHOOL 
MONTEREY.  CA  93940 


NAVAL  POSTGRADUATE  SCHOOL 

Monterey,  California 


THESIS 

HIGH-  RESOLUTION  COMPUTER  CALCULATION 

OF  OPTICAL  TRANSMITTANCE 

AT  SEA  LEVEL  OVER  MONTEREY 

by 

Nusret  Guner 

December   1978 

The 

sis 

Advisor:                   A.  W. 

Cooper 

Approved  for  public  release;  distribution  unlimited. 


T187436 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PACE  (Whmn  Dmtm  Enlmrmd) 


REPORT  DOCUMENTATION  PAGE 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


t.     REPORT  NUMBER 


2.  GOVT   ACCESSION   NO. 


3-      RECIPIENT'S  CATALOG   NUM8EB 


4.     TITLE   ;'mdiutim»j 


High  Resolution  Computer  Calculation 
of  Optical  Transmittance  at  Sea  Level 
Over  Monterey 


*•     TYPE  OF   REPORT   b  PERIOO  COVERED 

Master's  Thesis; 
December  1978 


«.  PERFORMING  ORG.  REPORT  NUMBER 


7.  AUTHORS 


Nusret  Guner 


•  •  CONTRACT  OR  GRANT  NL.M8ERC*; 


9.  PERFORMING  ORGANIZATION  NAME  AND  AOORESS 

Naval  Postgraduate  School 
Monterey,  California   93940 


10.     PROGRAM   ELEMENT.  PROJECT     TASK 
AREA  «   WORK  UNIT  NUMBERS 


II.     CONTROLLING  OFFICE  NAME   ANO  AOORESS 

Naval  Postgraduate  School 
Monterey,  California   93940 


12.  REPORT  OATE 

December  1978 


13.     NUMBER  OF   PAGES 
111 


14.     MONITORING  AGENCY   NAME  a    AOORESSf//  dIHsrwnt  Irom  Controlling  Olllco) 


IS.     SECURITY  CLASS,  (ot  thf  ripen) 

Unclassified 


1S«.     OECLASSIFI  CATION/' DOWNGRADING 
SCHEDULE 


14.     DISTRIBUTION   STATEMENT  (ot  thla  Rtport) 

Approved  for  public  release;  distribution  unlimited 


'7.     DISTRIBUTION  STATEMENT  (ot  tho  aostrmct  tntotod  In  Block  30,  tl  dUfmrmnt  troat  Rmpon) 


It,     SUPPLEMENTARY  NOTES 


19.     KEY  WORDS  (Contlnuo  on  r»r»»««  tido  U  r,»c»»tmrr  ond  Idontlfy  *y  block  numbor) 

Optical  Transmittance 
Molecular  Absorption 
Aerosol  Extinction 


20.     ABSTRACT  (Continuo  on  rovoroo  tido)  II  nocoooiy  •«*  idontitf  bf  block  ntmmkor) 

Computer  programs  have  been  developed  for  calculation  of  the 
nolecular  absorption  and  aerosol  extinction  for  horizontal  path 
propagation  of  radiation  in  the  wavelength  range  above  0.55  urn. 
The  AFCRL  line  parameters  compilation  is  used  as  input  data  with 
Local  weather  conditions  for  the  Monterey  Bay  and  model  of  typical 
aerosol  size  distribution  from  Shettle  and  Fenn.   Wind  speed  and 
relative  humidity  are  included  in  the  form  developed  by  Wells, 


DO 


FORM 
1  JAM  73 


1473 


EDITION  OP   1  MOV  *•  IS  OBSOLETE 
S/N    0102-014-4601   I 


UNCLASSIFIED 


SeCURlTY  CLASSIFICATION  OF  THIS  PAGE  (Whom  Dot*  Kntotod) 


UNCLASSIFIED 


(bcuwtv  ci a««i*ic ATiQM  or  tmis  »totrws«i  n*<«  r*«»»«<* 


ABSTRACT  (Cont'd) 

Gal  and  Munn.   Infinite  resolution  computations  have  been 
carried  out  for  1.06,  3.8  and  10.6  ym  and  results  tabulated 
for  monthly  and  yearly  average  weather  conditions  on  Monterey 
Bay.   Bandwidth- averaged  spectral  transmittances  have  been 
calculated  at  0.01  cm   and  0.5  cm  '  resolution  for  the 
same  conditions. 

Comparisons  with  previously  published  data  show  agreement 
to  1.8%  with  maritime  model  aerosol  Mie-scattering  calculations 
by  Selby  using  LOWTRAN  III  B  [Ref .  12] . 

Analysis  of  sensitivity  to  changes  of  pressure  and 
temperature  for  the  3.8  pm  range  showed  only  small  effects 
of  these  parameters.   However  a  change  in  water  vapor 
partial  pressure  results  in  a  0.75%/mbar  change  in 
transmittance  in  the  3.8  ym  region. 


DD      Form        1473  UNCLASSIFIED 

S/N    0102-014-6601  2        »icu«mtv  claudication  o'  tmii  p»ccr**««  o*f  *«r.».*> 


Approved  for  public  release;  distribution  unlimited. 

High  Resolution  Computer  Calculation 

of  Optical  Transmittance 

at  Sea  Level  Over  Monterey 

by 

Nusret  Guner 
Lieutenant,  Turkish  Navy 
B.S.,  Naval  Postgraduate  School,  1978 


Submitted  in "partial  fulfillment  of  the 
requirements  for  the  degree  of 


iMASTER  OF  SCIENCE  IN  PHYSICS 


from  the 

NAVAL  POSTGRADUATE  SCHOOL 
December  1978 


Q.I 


DUDLEY  KNOX  UBRARY 
NAVAL  POSTGRADUATE  SCHOOL 
MONTEREY.  CA  93940 


ABSTRACT 


Computer  programs  have  been  developed  for  calculation 
of  the  molecular  absorption  and  aerosol  extinction  for 
horizontal  path  propagation  of  radiation  in  the  wavelength 
range  above  0.55  urn.   The  AFCRL  line  parameters  compilation 
is  used  as  input  data  with  local  weather  conditions  for  the 
Monterey  Bay  and  model  of  typical  aerosol  size  distribution 
from  Shettle  and  Fenn.   Wind  speed  and  relative  humidity 
are  included  in  the  form  developed  by  Wells,  Gal  and  Munn. 
Infinite  resolution  computations  have  been  carried  out  for 
1.06,  3.8  and  10.6  ym  and  results  tabulated  for  monthly 
and  yearly  average  weather  conditions  on  Monterey  Bay. 
Bandwidth-averaged  spectral  transmittances  have  been  calcu- 
lated at  0.01  cm    and  0.5  cm    resolution  for  the  same 
conditions . 

Comparisons  with  previously  published  data  show  agree- 
ment to  1.8%  with  maritime  model  aerosol  Mie-scattering 
calculations  by  Selby  using  LOWTRAN  III  B  [Ref .  12] . 

Analysis  of  sensitivity  to  changes  of  pressure  and 
temperature  for  the  3.8  ym  range  showed  only  small  effects 
of  these  parameters.   However  a  change  in  water  vapor 
partial  pressure  results  in  a  0.75%/mbar  change  in  trans- 
mi  ttance  in  the  3.8  ym  region. 
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I.   INTRODUCTION 

A.   BACKGROUND 

The  light  bringing  us  information  from  a  distance, 
illuminating  our  surroundings  and  carrying  energy,  per- 
forms these  functions  only  by  travelling  in  the  atmos- 
phere.  Laser  light  in  the  form  of  a  highly  collimated 
beam  with  very  narrow  frequency  bandwidth  and  high  radiant 
power  is  equally  dependent  on  the  propagation  properties 
of  the  atmosphere.   The  characteristics  of  the  light: 
notably  its  color  (wavelength) ,  intensity,  polarization 
and  spatial  distribution  are  all  modified  by  interaction 
with  the  constituents  of  the  atmosphere  -  mostly  the  absorb- 
ing trace  gases  and  suspended  particles  in  the  size  range 
up  to  the  order  of  hundreds  of  micrometers. 

In  modern  military  technology  a  wide  range  of  systems 
are  in  use  or  under  development  for  surveillance,  imaging, 
target  detection,  designation,  identification  and  tracking, 
for  communications,  for  precision  guidance  of  munitions 
and  possibly  even  as  a  directed  energy  weapon.   These  oper- 
ate over  the  entire  "optical"  range  of  wavelengths,  from 
the  visible  or  near  ultra-violet  through  the  "far  infra  red" 
from  about  0.3  to  20  micrometers.   Operation  of  these  sys- 
tems is  limited  by  the  properties  of  the  atmosphere  as  a 
transmission  medium,  particularly  the  degradation  of 
resolution  by  scattering  from  particulates  and  turbulence, 
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and  loss  of  radiant  intensity  through  absorption  and  scatter- 
ing.  Due  to  strong  absorption  by  atmospheric  molecules, 
systems  are  constrained  to  operate  in  limited  wavelength 
ranges  or  "windows",  relatively  free  of  absorption  lines. 
Of  these  windows,  the  radiating  temperatures  of  many 
sources  of  military  interest  and  the  desirability  of  day/ 
night  operating  capability  have  concentrated  attention  on 
the  3-5  um  and  8-14  ym  infra-red  windows. 

In  particular  the  navy  has  a  special  interest  in  the 
DF  laser  as  a  candidate  moderate  or  high  power  laser  at  a 
wavelength  in  the  3-5  ym  range,  unaffected  by  the  high 
humidity  levels  found  in  the  naval  environment.   Of  particu- 
lar interest  is  the  3.8007  P~(8)  transition  line  which 
shows  up  to  9  6%  transmission  over  a  kilometer  path  at  sea 
level. 

The  need  to  predict  the  operational  performance  of 
these  systems  leads  to  a  requirement  to  predict  the  spec- 
tral transmittance  of  the  atmosphere  in  the  infra-red,  as 
a  function  of  the  meteorological  parameters  which  determine 
the  absorption  and  scattering  levels  —  the  pressure,  tem- 
perature, air  composition,  vertical  gradients,  number  den- 
sity composition  and  size  distribution  of  suspended  particles 

Large  scale  measurement  programs  have  been  undertaken 
in  many  countries  with  the  objective  of  building  a  large 
library  of  meteorological  and  transmission  data.   Consider- 
able effort  (notably  at  AFGL)  has  been  devoted  to  the 
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development  of  computer  programs  for  transmittance  calcu- 
lation based  on  weather  data,  and  the  establishment  of  a 
number  of  "standard"  atmospheric  models. 

The  current  modelling  effort  was  undertaken  in  support 
of  the  experimental  program  of  transmittance  measurement 
on  ranges  over  Monterey  Bay  being  carried  out  by  the 
Environmental  Physics  Group  at  the  Naval  Postgraduate 
School.   This  work  involves  measurement  in  all  the  suitable 
atmospheric  windows  from  the  visible  to  14  urn,  with  both 
laser  and  broadband  sources.   The  calculational  program 
has  therefore  been  directed  to  computation  of  molecular 
absorption  and  aerosol  extinction  coefficients  for  sea 
level  propagation  under  a  range  of  weather  conditions  sta- 
tistically typical  of  Monterey.   This  has  been  done  for  a 
variety  of  wavelengths  appropriate  to  the  measurements , 
and  at  high  resolution  for  use  with  laser  sources,  and 
averaged  over  filter  bandwidths  for  broadband  sources.   The 
Air  Force  Cambridge  Research  Laboratory  Atmospheric  Absorp- 
tion Line  Parameters  Compilation  [Ref.  7]  in  magnetic  tape 
form  is  used  as  input  data  for  the  absorption  calculations. 
A  Mie  scattering  program  based  on  the  Van  De  Hulst  approxi- 
mation with  the  model  particle  distribution  of  Shettle  and 
Fenn  [Ref.  8]  has  been  used  for  the  aerosol  scattering. 
Wind  speed  and  relative  humidity  effects  have  been  included 
using  the  drop  size  growth  model  of  Wells,  Gal  and  Munn 
[Ref.  18].   A  sensitivity  analysis  has  been  carried  out  with 
respect  to  pressure  and  temperature. 
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Basically,  this  thesis  describes  a  program  to 

a)  Develop  a  computer  code  to  calculate  atmospheric 
transmittance  for  specific  wavelengths  in  the  near  and 
mid  IR. 

b)  Develop  a  code  to  give  transmittance  and  molecular 
absorptance  with  bandwidth  appropriate  to  grey-body  sources 
in  the  near  and  mid  IR. 

c)  Apply  these  to  prediction  of  transmittance  over 
marine  optical  paths,  over  Monterey  Bay,  based  on  the 
available  meteorological  data  base  for  the  area,  for  con- 
ditions appropriate  to  the  optical  propagation  experimental 
measurement  program. 

Curves  have  been  prepared  giving  predicted  molecular 
and  aerosol  extinctions  for  statistical  monthly  and  yearly 
average  meteorological  conditions  on  Monterey  Bay. 

Further  sections  of  this  thesis  describe  the  theoreti- 
cal background  and  the  details  of  the  algorithms,  and  provide 
tabular  and  graphical  results  of  the  calculations.   Computer 
program  listings  are  given  in  the  appendices. 

B.   ATMOSPHERIC  PHENOMENA 

Optical  transmission  of  energy  is  usually  described 
by  the  transmittance, 


t   =   exp(-yL)   =   ~[^j-  (1-1) 
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where  y  is  called  the  extinction  coefficient,  and  L  repre- 
sents the  length  of  the  path  which  the  light  has  traveled. 
This  expression  follows  from  Beer's  Law,  which  for  linear 
propagation  of  monochromatic  radiation  at  frequency  v  has 

the  form 


37  I(z,v)   =   -y  (v)  I(z,y)  (1-2) 


Here  I(z,y)  represents  the  intensity  of  monochromatic 
radiation  at  a  distance  z. 

Beer's  Law  in  differential  form  (1-2)  may  also  be 
written  as 


I(L,y)   =   1(0, y)  e"u(v)L  (1-3) 


for  a  given  path  length  L.   The  transmittance  t  over  a 
path  of  length  L  is  then  given  from  the  equation  (1-3) 


T(v)   =  Kiwv)   =  e-y(v)L 

T^v;      l(0,v)      e 


which  is  the  same  as  equation  (1-1)  for  monochromatic 
radiation. 

The  extinction  coefficient  y  is  generally  the  result 
of  two  processes:   absorption  and  scattering.   Both  mole- 
cules and  particles  suspended  in  the  atmosphere  cause 
absorption  and  scattering.   It  may  be  useful  to  define  the 
extinction  coefficient  y  as  a  sum  of  four  terms 
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V     =  km  +  2m  +  k  +  Sa  (1-4 

m    m    a    a 


where  k  ,  3  ,  k  ,  3^  are  molecular  absorption,  molecular 
m   m   a    a  c 

scattering,  aerosol  absorption,  aerosol  scattering 

coefficients . 

An  aerosol  extinction  coefficient  u   is  sometimes  used 

a 

for  the  summation  of  aerosol  absorption  and  scattering. 

Similarly  a  molecular  extinction  coefficient  u   is  used 

m 

for  the  summation  of  molecular  absorption  and  scattering. 

From  equations  (1-1)  and  (1-4)  the  total  transmittance 
through  a  path  in  the  atmosphere  may  be  written  as 


i (total)   =   t (molecular  absorption)  * 
x (molecular  scattering)  x 
t  (aerosol  absorption)  * 
t (aerosol  scattering) 

where  x (molecular  absorption)  may  be  written  in  terms  of 
the  atmospheric  components  as 

t (molecular  absorption)   =   t (uniformly  mixed  gases  line  abs.)  x 

x  (ozone  line  abs.)  * 
x (water  vapor  line  abs.)  x 
t (water  vapor  continuum  abs . )  * 
x (nitrogen  continuum  abs.) 
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Although  the  molecular  scattering,  aerosol  absorption 
and  aerosol  scattering  coefficients  are  smooth  functions 
of  frequency  of  monochromatic  radiation,  the  molecular 
absorption  coefficient  is  strongly  dependent  on  frequency. 

Molecules  which  cause  absorption  and  scattering  of  a 
light  beam  in  the  atmosphere  may  be  divided  into  two  types: 
uniformly  and  non-uniformly  mixed  gases.   Carbon  dioxide, 
nitrous  dioxide,  carbon  monoxide,  methane,  oxygen  and 
nitrogen  are  assumed  as  uniformly  mixed  gases  for  which 
their  mixing  ratios  by  volume  with  each  other  are  constant, 
independent  of  pressure  and  temperature.   These  mixing  ratios 
are  tabulated  in  table  4.   Water  vapor  and  ozone  are  con- 
sidered as  non-uniformly  mixed  gases,  whose  relative  con- 
centrations must  be  specified. 

Six  different  models  of  the  atmosphere  are  frequently 
used  in  calculations;  these  are  the  tropical,  midlatitude 
summer,  midlatitude  winter,  subarctic  summer,  subarctic 
winter  and  U.S.  Standard  atmosphere  196  2  model  described  by 
McClatchey  and  D ' agati  [Ref.  4].   Pressure,  temperature, 
density,  water  vapor  and  ozone  concentrations  in  the  air 
for  each  model  at  sea  level  are  tabulated  in  table  5 . 

The  number  densities  of  aerosol  particles  which  also 
cause  absorption  and  scattering  of  a  light  beam  are  more 
complicated  to  determine.   It  is  necessary  to  know  the  size 
distribution  of  aerosol  particles  in  addition  to  the  number 
density,  in  order  to  compute  the  aerosol  scattering  and 
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absorption  coefficients.   The  size  distribution  function 
is  denoted  by  n(r)  and  given  as 


»<r)   -  ^l£l  (1-5, 


where  N(r)  is  called  the  number  density  which  is  the  number 
of  aerosol  particles  with  radius  less  than  r  per  unit  volume 

Several  different  analytical  models  have  been  developed 
for  the  size  distribution.   A  "continental"  aerosol  model, 
which  includes  both  rural  and  urban  aerosol  conditions, 
and  a  "maritime"  aerosol  model  are  used  for  the  lower 
atmosphere. 

The  rural  aerosol  model  used  by  Shettle  and  Fenn  is 
assumed  to  be  composed  of  a  mixture  of  70%  water  soluble 
substance,  consisting  of  ammonium  and  calcium  sulfate  and 
organic  compounds,  and  30%  dustlike  aerosols  [Ref.  8]. 

In  urban  areas  the  rural  aerosol  background  is  modified 
by  the  addition  of  aerosols  from  combustion  products  and 
industrial  sources.   The  proportions  of  soot-like  aerosols 
and  the  rural  aerosol  type  aerosol  mixture  are  assumed  to 
be  35%  and  65%  respectively  [Ref.  8] . 

The  maritime  aerosol  model  represents  the  aerosol  com- 
position and  size  and  number  distributions  over  the  oceans. 
It  is  significantly  different  from  the  continental  aerosol 
types.   The  maritime  model  is  due  to  salt  particles  caused 
by  the  evaporation  of  the  seaspray  droplets  also  added  to 
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the  rural  type  aerosols,  from  which  most  of  the  very  large 
particles  have  been  eliminated. 

The  general  form  of  the  size  distribution  for  these 
models  is  given  by  Shettle  and  Fenn  [Ref.  8]  as 

2        N.  (log  r  -  log  r. ) 2 

n(r)   =    I    ( ± T7T)exp[ * 1 — ] 

i=1  (In  10)rai(2Tr)x/z  2c/ 

(l-6a) 

where  the  parameters  in  the  equation  above  are  tabulated 
in  table  6. 

The  size  distribution  given  by  equation  (l-6a)  is 
relative  humidity  and  wind  speed  independent.   Since  the 
number  of  oceanic  particles  with  greater  radius  in  the  size 
distribution  increases  with  increasing  relative  humidity  and 
wind  speed,  a  modified  maritime  model  has  been  proposed 
for  calculation  of  aerosol  extinction  coefficients  taking 
care  of  the  growth  of  particle  radius  [Ref.  18].   This 
gives 


n(r)   =   p-{0.47(^r)  4  +  2  .  3MRa  (~)  exp  [-8  .  5  (p-)  Y]  }   (l-6b 


where  C  is  a  normalization  constant,  and  M   is  the  mixing 
ratio  of  seaspray  aerosols  with  the  rural  aerosols.   Or 
more  clearly,  the  mixing  ratio  of  rural-oceanic  aerosols 
in  the  equation  (l-6b)  is  1:MR- 

Particle  growth  is  included  through  the  growth  factor  F 
which  adds  the  effect  of  relative  humidity  into  the  equation 


(l-6b) .   It  is  given  as 


F   =   1-0.9  In (1- 


100 


where  R.H.  is  the  percent  relative  humidity.   The  effect 
of  F  is  to  increase  the  drop  radius  to 


r      =   F  r 
R.H.         o 


where  r   is  the  radius  for  zero  percent  relative  humidity. 

O  sr  j. 

F  is  normalized  to  80%  relative  humidity,  given  F'.   There- 
fore factor  F'  in  the  equation  (l-6b)  is 


,      F(R.H.) 
F(80) 


Y  in  equation  (l-6b)  describes  the  dependence  on  the  wind 
speed.   It  has  the  form 


Y   =   0.384  -  0.00293  V1'2^ 


where  V  is  the  wind  speed  in  meters  per  second.   a  in 
equation  (l-6b)  also  represents  the  wind  speed  dependence 
and  it  is  given  as 


a   =   250  +  750  V1"16  ,    when   V  <  7  m/sec. 


a   =   6900  v      ,         when   V  >  7  m/sec 
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So  the  inputs  for  equation  (l-6b)  are  wind  speed  V  and 
relative  humidity  R.H.   The  normalization  constant  C  may 
be  determined  from  the  relation  between  visibility  and 
aerosol  extinction  coefficient  at  0.55  micrometers. 

C.   HF/DF  LASER  AT  THE  NAVAL  POSTGRADUATE  SCHOOL 

Although  this  thesis  will  deal  with  propagation  at 
sea  level  at  any  frequency  of  radiation,  primary  interest 
will  be  the  frequencies  of  the  HF/DF  laser.   These  fre- 
quencies are  tabulated  in  tables  7  and  8 . 

To  be  used  for  the  propagation  measurement  program  at 
the  Naval  Postgraduate  School,  a  small  scale  hydrogen 
fluoride  (HF)  or  deuterium  fluoride  (DF)  laser  was  located 
on  the  roof  of  Spanagel  Hall.   The  system  consists  of  a 
vacuum  system,  a  power  supply  and  the  laser  head.   The 
laser  is  centered  around  a  laser  head  which  contains  dis- 
charge tube,  mixing  chamber  or  cavity  and  resonant  cavity. 
The  system  uses  four  gases:   a)  Sulfur  hexaf luoride,  b)  oxy- 
gen, c)  helium,  d)  hydrogen  or  deuterium.   A  mixture  of 
sulfur  hexaf luoride ,  oxygen  and  helium  is  injected  into  the 
discharge  tube  at  the  end  opposite  the  mixing  chamber. 
Then  the  sulfur  hexafluoride  is  dissociated  by  electron 
impact  in  the  discharge  tube  by  making  use  of  an  electric 
discharge.   Oxygen  helps  the  sulfur  hexafluoride  to  dis- 
sociate into  free  fluorine  radicals  by  the  reaction 


SF„  +  0-   -   S0„Fo  +  2F 
4     2        2  2 
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The  hydrogen  or  deuterium  is  injected  perpendicular  to  the 
flow  in  the  mixing  chamber.   Then  chemical  reaction  takes 
place  in  this  mixing  region  forming  the  vibrationally 
excited  HF  or  DF  molecules. 

The  optical  resonator  of  the  system  consists  of  a 
pair  of  Brewster  windows  and  mirrors.   The  Brewster  windows 
which  linearly  polarize  the  output  are  chosen  from  two  inter- 
changeable pairs,  one  of  BaF~  (0.25-15  ym  transmitting) 
and  the  other  sapphire  (0.14-6.5  \im   transmitting). 

The  system  has  two  options  for  operation:   a)  Multi- 
line operation,  b)  Single  line  operation.   The  output 
mirror  should  be  replaced  by  the  diffraction  grating  for 
single  line  operation,  while  two  mirrors  (one  partially 
transmitting)  are  used  for  multiline  operation. 

Optimized  power  output  from  lasing  of  hydrogen 
fluoride  was  reported  [Ref.  9]  as  6.5  Watts  with  the  mass 
flow  rates 


SF    =   0.600  grams/sec. 
6 


O-p    =   0.182  grams/sec 


B.  =   0.019  grams/sec. 

He    =   0.053  grams/sec. 


Optimized  multiline  power  output  of  deuterium  fluoride 
lasing  was  also  reported  [Ref.  9]  as  2.5  Watts  with  the 
mass  flow  rates 
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SF6  =  0.680  grams/sec. 

02  =  0.270  grams/sec. 

D~  =  0.043  grams/sec. 

He  =  0.063  grams/sec. 

In  the  case  of  single  line  operation,  the  wavelengths 
and  corresponding  relative  intensities  of  fourteen  individual 
DF  and  fifteen  individual  HF  lasing  transitions  are  shown 
in  figure  3  and  figure  4. 

The  block  diagram  of  the  main  elements  of  the  system 
is  also  shown  as  figure  5 . 
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II.   OPTICAL  PROPAGATION  THROUGH  THE  ATMOSPHERE 

A.   MOLECULAR  SCATTERING  COEFFICIENT 

The  simple  Lorentz-type  molecular  model  presented  here 
is  shown  in  figure  1.  It  is  assumed  nonionized,  nonpolar, 
isotropic,  linear  and  lightly  damped. 

The  model  contains  molecular  mass  residing  in  the 
atomic  nuclei  and  carrying  a  net  positive  charge.   Orbital 
electrons  shared  by  the  constituent  atoms  are  represented 
by  a  concentric  shell,  providing  a  balancing  negative 
charge.   Force  directed  toward  the  central  mass  is  char- 
acterized by  a  spring  constant  k.   Let  m  be  the  electron 
mass.   If  the  system  is  disturbed  by  an  electric  field  E 
in  the  z  direction,  then  the  equation  of  motion  may  be 
written  for  a  displaced  electron  as 


d2z 
-  eE   -  kz   =   m  — *•  (2-1) 

dt2 


where 


E   =  E    sin  cot. 
z    oz 


Equation  (2-1)  can  also  be  written  as 


2  -eE 

14  +   <A      =      —  sin   cot  (2-2) 

dt2 
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1/2 
where  uq   =  (k/m)     is  defined  as  the  resonance  frequency, 

corresponding  to  a  molecular  transition. 

The  steady  state  solution  to  equation  (2-2)  is 

_       ~e      t-, 

z   -   s ^—  E    sin  cot 

z  ^     A   oz 
m(co   -  co  ) 

Now,  the  oscillator  dipole  moment  can  be  written 
as 

2  , 
„             e  /m 
p   =   -ez   =  —*£ s-  E    sin  oot 

2     2   oz 

CO    -  CO 

o 

If  we  denote  the  maximum  value  of  the  dipole  moment 

as  p  ,  then  it  is 
o 

g  /m  n  r 

ro       2     2   oz 

co   -  CO 

o 

The  oscillating  dipole  moment  has  maximum  value  p 
given  by  equation  (2-3)  and  produces  secondary  waves. 
Reradiation  from  the  oscillating  dipole  then  produces 
electric  and  magnetic  fields  at  every  point  in  the  coor- 
dinate system  shown  in  figure  2.   In  the  short  dipole 
approximation  (i.e.,  radius  of  scatterer  is  much  less  than 
wavelength  of  incident  wave)  and  far  field  (i.e.,  distance 
to  scatterer  is  much  larger  than  radius  of  scatterer)  the 
electric  and  magnetic  field  vectors  can  be  written  as 
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2 

-*■  to   P  • 

£     -  o  ysin    9       .     ..         ^      ~ 

E      -        4tt  — sm(kr-oot)    9 


2 

-*■             to    p  . 

□      -            o  ysin   9           ,.         .. 

B     ~     TT  vr        sin(kr-o)t) 


where 


v   is  the  velocity  of  wave, 

y   is  the  magnetic  permeability, 

9 ,  (p   and  r  are  the  variables  in  spherical  coordinates 

The  poynting  vector  denoted  by  S  describes  scattered  power 
per  unit  area,  and  is  given  here  as 

+   +  2 

J     EXB      y,w  Po,2.sin  9,2  .  2 ,,      *    ,0  . , 
S   =  -y-     =     -{—££-)     ( — - — )  sin  (kr-ut)r    (2-4) 

Scattered  power  per  unit  solid  angle  can  be  found  from 

equation  (2-4),  first  taking  the  time  average,  then  multi- 

2  . 

plying  by  r  .   This  gives  the  scattered  intensity  as 

2 

I(9)       "      2^(_4l~)     °in    (8)       "      2^{~~~2 2~]    Eoz 

4  7Tm(cu   -  oj  ) 

(2-5) 
the  power  of  the  incident  wave  contained  in  unit  area  can 
be  given  as 


<S>.   •  ,   .   =   i  v  e  E*  (2-6) 

incident      2       oz 


where  z    is  the  permittivity  of  the  medium. 
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The  angular  cross  section  of  a  molecule,  denoted  by 
am(9),  may  be  defined  as  that  cross  section  of  incident 
wave  acted  on  by  the  molecule,  having  an  area  such  that 
power  flowing  across  it  is  equal  to  the  power  scattered 
by  the  molecule  per  unit  solid  angle  at  an  angle  9.  It 
is 


I     ( fi ) 
a    /a)       -        scat _   equation  (2-5) 

<S> incident     equation  (2-6) 


am(9)   ~   {—]     {  —  2 27}    "  {-4T]     (E~ }       (2"7) 

m (oj   -  a)  )  oz 


It  is  necessary  to  define  the  polarization  P  in  order  to 
find  a  suitable  form  for  the  second  term  in  the  right  hand 
side  of  equation  (2-7).   P  is  defined  as  total  dipole 
moment  per  unit  volume  and  it  is  given  as 

P   =  Np   =  ve  E  «  =   (K-l)e  E  .. 

^o      A  o  eff  o  eff 

for  a  linear  and  isotropic  medium,  where 

N     is  number  of  molecules  per  unit  volume, 

X     is  electric  susceptibility, 

K      is  dielectric  constant, 

E  ,.*■      is  the  effective  electric  field  which  is 
eff 

E  +    p 


oz    3e 
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Now  we  can  write  P  in  terms  of  K  and  E   ,  as 

oz 


3  K-l  £ 

P   =   2.  E 

K+2      oz  ' 


from  which 


PQ      3(X-l)e 


E         N(K+2) 
oz 


The  index  of  refraction  is  given  here  as 


n   =   7:   =   /KK    =   /K 
v        m 


since  K   =  1  for  a  nonmagnetic  medium. 
m  3 

Now  equation  (2-7)  can  be  written  as 


2   3e  (n2-l)  - 

a(8)   =   (J^L-)(_°  )2 

47T    N(n2+2) 


a  ,  defined  bv  the  equation  below,  is  called  the  total 
m 

scattering  cross  section. 


4tt 

a    =    (      a    (Q)    dQ. 
m      ;    m 


or 


a   =  -il4(3(f-1})2  (2-8 

m      3N  r  (n  +2) 
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where 

2ttv 


CO    = 


A 

was  substituted. 

The  total  scattering  cross  section  given  by  equation 
(2-8)  is  defined  similarly  to  the  angular  scattering  cross 
section,  as  the  power  scattered  in  all  directions. 

For  air  the  refractive  index  is  very  close  to  one. 

2 

Therefore  n  +2  is  nearly  equal  to  three,  and  equation  (2-8 

may  be  written  as 


8  tt   ,  2  .  .  2  ,  _  _ 

Qm   =   — T-?(n  ~1)  (2'9 


This  form  of  the  total  scattering  cross  section  (2-9) 
has  been  derived  for  a  single  molecule  disturbed  by  incoming 
electromagnetic  radiation.   It  has  the  units  of  area.   For 
many  molecules  in  a  unit  volume,  the  total  scattering 
coefficient  3   may  be  defined,  assuming  each  molecule  has 
the  same  total  scattering  cross  section, 


3 

!    =   Na    =   -^-(n2-!)  (2-10 

m        m      3NA4 


The  molecular  scattering  coefficient  (2-10)  has  the  units 
of  1/length. 

Since  atmospheric  gases  are  not  quite  isotropic,  a 
depolarization  factor  P,  is  introduced  [Ref.  1]  and  inserted 
into  equation  (2-10) .   The  molecular  scattering  coefficient 
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may  then  be  written  as 


«tt3   9    9  6  +  3P,   _  T 


taking  into  account  the  pressure  and  temperature  dependence. 
Here  PQ  and  TQ  refer  to  standard  conditions,  1013  mb , 
273  K. 

Assuming  air  to  be  an  ideal  gas  and  using  the  ideal 
gas  law,  air  is  found  to  have  a  mole  density  of  44.631  moles 
per  cubic  meter  at  a  pressure  of  1013  millibars  and  tempera- 
ture of  273  degrees  Kelvin.   This  is  equivalent  to 

25 
2.688  *  10    molecules  per  cubic  meter. 

For  a  depolarization  factor,  a  value  of  0.035  is 

used,  following  Penndorf  [Ref.  1] .   Although  the  index  of 

refraction  of  air  depends  on  the  wavelength,  this  variation 

2    2  -7 

is  small.   Therefore  (n  -1)   is  taken  as  3.0884  x  10 

which  corresponds  to  a  wavelength  0.55  micrometers  [Ref.  2]. 
Using  these  parameter  values  given  above,  the  atmospheric 
molecular  scattering  coefficient  may  be  written  as 


3m  -  10-°y  10'20.^3)  <^»  -its  of  !/*» 


where 

a  has  units  of  cm. 

P  has  units  of  mbars . 

T  has  units  of  degrees  Kelvin 
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McClatchey  and  d'Agati  [Ref.  4]  have  used  the  slightly 
different  form  of  the  molecular  scattering  coefficient, 


_   9.807  x  iq"20   P    273,  „  .., 

Sm  ~   A4.0117 (I0T3)  (~ }  (2_11) 


which  was  selected  as  giving  a  "best  fit"  to  the  experi- 
mental data  of  Penndorf  [Ref.  6], 

B.   MOLECULAR  ABSORPTION  COEFFICIENT 

The  model  to  derive  the  molecular  absorption  coeffi- 
cient is  the  same  as  the  one  we  developed  in  the  previous 
section,  but  now  it  is  heavily  damped.   When  the  model  is 
disturbed  by  an  electric  field,  the  electron  moves  with  a 
damping  coefficient    and  restoring  force  strength  k.   The 
equation  of  motion  may  be  written  in  differential  form  such 
that 


d2r      dr   ,  +  t  £,*  -iwt 

m  — j  +  my   ^—  +  kr   =   -eE   =   -eE(r)e 

dt       at 


The  steady  state  solution  to  this  equation  is 

-e/m   Z, 


r  = 


2     2 
0)   —  0)   —  iy^> 

o 


The  polarization  due  to  N  molecules  per  unit  volume  is 


-r 


P   -   -eN:  -  ^ lT—   E.ff        (2-12) 
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where  we  again  replaced  the  applied  electric  field  with 
an  effective  electric  field  to  account  for  polarizability 
of  the  medium. 

After  substituting  Eeff  =  E  +  P/3eo,  the  polarization 
becomes 


p   -    e  N/m 


2     2 

co,  -  co   -  iyco 


where 


2     2     e2N 

CO    =  CO    - 

1     o    3£  m 
o 


But  the  polarization  is  also  given  as 


2 


P  -   xeo  Eeff  =  (K-l)  eQ  Eeff  =  (n  -1)  eQ   Eeff 

(2-13) 
Combining  equation  (2-12)  and  (2-13)  one  can  get 

2  , 

P   1/2  e  N/£om     1/2 

«      =      (1+^E}     =  {1+- T— } 

o  co,  -  co   -  ly  cu 

(2-14) 

for  the  index  of  refraction. 

A  more  general  form  for  equation  (2-14)  is  written 
below  with  j  referring  to  different  types  of  molecules 
with  oscillator  strength  f ■ . 


o2       Nifi         i  /-> 

,-,  e      V  ^    J  x  1/2  ,~    ,  -  X 

n    =    (1  +  nz  l  — 2 — — }  (2~lD) 

o  co,   -  co   -  iyco 
j 
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Near  a  resonance  frequency,  equation  (2-15)  can  be 
approximated  by 


2  2    2  ? 

eHS.f  .     (a?,  -co  )         ezN  f 

n   =   [l  +  1-1  13 ]  +  r    j  3   ^1 


2s:  m     .2    2,2^  2  2J   L2e  m      2    2,2^  2  2 

O      (go,  .-co  )  +y  w        O      oo,-j-oo  )  +y  oj 

(2-16) 


using  the  Binomial  Expansion,  since  n  is  very  close  to 
unity. 

A  plane  wave  going  in  the  z  direction  with  wave  vector 
K  may  be  written  as 


E   =   E   ei(Kz-aJt)  (2-17) 

o 


The  relation  between  wave  vector  and  refractive  index 
is 


K   =  -  £  n  (2-18) 

c 


K  is  complex,  since  n  is  complex,  so  defining  K  and  n 

K  =   3  +  ia  (2-19) 

n   =   n  +  ik  (2-20) 

where  S,  a,  n,  k  are  the  so-called  propagation  constant, 
absorption  coefficient,  index  of  refraction,  and  extinction 
coefficient  respectively. 
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Using  equations  (2-16,18,19,20),  it  can  be  shown  that 

e2N,f .         2  2 

2e  mc  y  7^ 2,2^  2  2  (2  21) 

X    (oo,  .-co  )  +co  y 


Substituting  K  from  equation  (2-19)  into  equation  (2-17) 
one  can  get 


£     =   E   e"az  ei(3z"ajt) 
o 


The  electric  field  decreases  by  a  factor  of  e~az, 
and  the  intensity  falls  off  by  a  factor  of  e     along  the 
path  length  z.   Since  a  corresponds  to  absorption,  many 
workers  use  the  so-called  Lorentzian  line  shape  given  in 
equation  (2-21)  in  order  to  find  an  expression  for  molecular 
absorption  coefficient. 

Substituting  the  expressions  for  frequency  oj  =  2ttv    and 
oj,  ■  =  2pv   into  equation  (2-21)  a  becomes 


a   = 


[ — (v  -v)  ]   +  1 
Y   o 


or 


a   = 


A  Av2 


2      2 

(v  -v)   +  Av 


(2-22 
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where  A  is  a  constant,  Av  is  y/4tt  and  defined  as  half 
width  at  half  maximum. 

If  a  is  normalized  in  the  range  from  minus  infinity  to 
plus  infinity,  the  expression  for  the  Lorentzian  line 
shape  becomes 


g(v,   =   Av 


2      2    ' 


describing  the  spectral  distribution  of  absorption. 

Finally,  the  molecular  absorption  coefficient  can  be 
written  as 


K  (v)   = 


m  r  ,     .  2  ,  27    ' 

it  [  (v  -v)   +  a  ] 

provided  S  =  /k  (v)dv  for  single  absorbing  line,  where 

S   is  line  intensity  per  absorbing  molecule  in 
the  units  of  cm"v  (molecules/cm^)  , 

a   is  Av  half  width  at  half  maximum  in  the  units 
of  cm"S 

v   is  incoming  wave  frequency  in  the  units  of  cm   , 

v    is  resonance  frequency  for  absorbing  line  in 
the  units  of  cm 


The  radiant  frequency  v  is  unique:   this  gives  absorption 
as  a  spectral  quantity.   Radiation  from  laser  sources  has 
very  narrow  line  width,  which  may  be  assumed  infinitely 
narrow  for  calculation  purposes.   The  resonance  frequency 
for  each  absorbing  line  is  independent  of  pressure  and 
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temperature.   The  line  width  is  assumed  to  be  a  pressure 
broadened  Lorentzian,  with  its  pressure  and  temperature 
dependence  given  by  the  expression  [Ref.  7] 

a   =   <x(PQ/Ts)  f-    (^-)1/2  (2-23) 

o 


where 


P   =   1  atm.  and 

T    =   296  degrees  Kelvin. 

The  total  intensity  of  any  line  is  pressure  independent 
and  its  temperature  dependence  is  given  by  the  expression 
[Ref.  7] 


Q  (T  )   Q  (T  )  T-T 

S(T)   =   S(Tg)  q  (T)   xq  (T)    xexp[1.439E"(-^-)]    (2-24) 


where  QTT  and  Q   are  vibrational  and  rotational  partition 
functions  respectively,  E"  is  the  lower  state  energy  of 
the  transition  corresponding  to  that  line,  in  units  of  cm 

Rotational  partition  functions  are  also  temperature 
dependent,  of  the  form 


W       ,T,j 

Qr(T)        4sj 


The  corresponding  values  of  Qv  and  j  are  tabulated  for  each 
molecule  in  the  air  [Ref.  7]  and  presented  in  table  2. 
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Finally  the  total  molecular  absorption  coefficient  due 
to  water  vapor,  ozone,  carbon  dioxide,  methane,  carbon 
monoxide,  oxygen,  and  nitrous  oxide  in  the  air  is  given 
by 


S  .  . a .  .m . 
Km   =  H       r^  XM 2- 


j  i  ^^ij)   +  *   1 


13       i]' 

Summation  over  j  corresponds  to  different  types  of  mole- 
cules, and  summation  over  i  corresponds  to  different  lines 
for  a  given  molecule. 

Values  of  v.  ,  a   (P  ,T  ) ,  S..(T  )  have  been  compiled 
for  100,000  lines  between  1  ym  and  the  far  infrared  by 
McClatchey  and  others  at  Air  Force  Cambridge  Research 
Laboratories  and  are  available  on  magnetic  tape  [Ref.  7]. 
Also  a  new  version  of  the  compilation  including  over 
139,000  lines  between  0.2  ym  and  30  urn  was  recently  reported 
by  Air  Force  Geophysics  Laboratory  [Ref.  3]. 

Calculation  for  K  at  any  pressure  and  temperature  may 
be  made  by  using  equations  (2-23,24,25)  with  the  data  from 
this  tape. 

m.  in  equation  (2-25)  represents  the  number  of  j  type 

molecules  per  unit  volume.   The  resultant  K   is  in  the  units 
c  m 

of  1/length  and  generally  1/kilometer  is  used. 

It  should  be  noted  that  the  Lorentzian  line  shape  must 
be  modified  for  pressures  less  than  100  millibars;  at  pressures 
between  10  and  100  millibars  doppler  broadening  becomes 
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dominant  [Ref.  7],  while  at  pressures  less  than  10  milli- 
bars, the  Voigt  profile  should  be  used.   The  Lorentzian 
line  shape  is  good  enough  to  be  used  at  low  altitudes,  but 
line  wings  should  be  truncated  at  20  cm   from  the  line 
centers,  since  uncertainties  begin  to  arise  beyond  that 
range. 

A  computer  program  calculating  the  molecular  absorption 
coefficient  is  prepared  and  presented  as  Appendix  A.   This 
program  also  includes  the  water  vapor  continuum  absorption 
calculation  which  will  be  explained  in  the  next  section. 
The  algorithm  for  the  program  is  also  discussed  in  Section 
III. 

C.   CONTINUUM  ABSORPTION  COEFFICIENT 

Continuous  regions  of  absorption  in  the  atmospheric 
windows  occur  near  250  0  cm    (4  ]_im)  and  from  7  00  cm    to 
1250  cm"1  (8-14  ym) . 

In  other  spectral  regions ,  the  contribution  from  nearby 
absorption  lines  is  much  greater  than  that  from  the  continuum 
absorption,  so  that  for  practical  purposes  the  continuum 
effect  can  be  neglected,  although  it  may  be  greater  than 
in  the  windows. 

Extreme  wings  of  strong  collision  broadened  absorption 
-lines  centered  more  than  20  cm   away  or  pressure  induced 
absorption  resulting  from  transitions  that  are  forbidden  for 
unperturbed  molecules  may  cause  continuous  absorption.   In 
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the  8-14  um  region,  existence  of  the  water  dimer  (H  0:H.O) 
may  also  cause  continuous  absorption. 
1.   Water  Vapor  Continuum 

The  absorption  coefficient  due  to  the  water  vapor 
continuum  can  be  described  empirically  by  the  expression 
[Ref.  4] 


Km        =   [Cs(v,T)Ps  +  C  (v,T)P  ]W  (2-26 

cont . 


where  C_  is  a  self  broadening  coefficient  due  to  self 
broadening  from  collisions  of  water  molecules  with  other 
water  molecules,  CN  is  a  foreign  gas  broadening  coefficient 
due  to  nitrogen  broadening  from  collisions  of  water  mole- 
cules with  air  molecules  which  are  mostly  nitrogen. 

C   and  C   are  both  temperature  and  frequency  dependent, 
a.   8-14  ym  region  [Ref.  4] 

C   is  given  by  the  expression 


CN(v,T)   =   0.002  x  Cs(v,296 


C   is  given  by  the  expression 
S 


Cs(v,T)       =      Cs(v,296)     x    exp[1800(i  -   j~0 


where    Cs (  v  ,296)     is    given   as 
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C  (v,296)   =   4.18  +  5578  x  exp(-7.87  x  lo"J  v) 


where  Cg  is  in  the  units  of  (precipitable-cm) ~1  atm"1  or 
(pr.  cm)    atm   which  is  (grams/cm  )    atm   . 
b.   3.5-4.2  ym  region  [Ref.  4] 

Cs  and  C^  are  given  by  the  expressions  below; 


CN(v,T)   =   0.12Cs(v,T) 


Cs(v,T)   =   Cs(v,296)  x  exp[1350(^  -  ~g)] 


where  Cs  has  the  units  of  (pr.  cm)    atm   .   The  values  of 

Cs(v,296)  at  different  frequencies  are  reproduced  [Ref.  4] 

and  tabulated  in  table  3.   W  in  the  equation  (2-26)  has 

the  units  of  (pr.  cm)  per  kilometer  which  is  the  same  as 

grams  per  square  centimer  per  kilometer.   Then  the  water 

vapor  continuum  absorption  coefficient  K       can  be 

_-,  cont. 

obtained  in  (kilometers) 

2 .   Nitrogen  Continuum 

Although  the  nitrogen  continuum  is  not  included  in 
this  thesis,  the  magnitude  of  this  effect  is  of  the  order 
of  0.001  km   at  a  wavelength  of  3.8  urn.   The  maximum 
effect  which  is  of  the  order  of  0.1  km    occurs  at  a  wave- 
length of  4.3  urn.   This  effect  may  be  neglected  outside 
the  3.7-4.7  ym  region  [Ref.  4]. 
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D.   AEROSOL  EXTINCTION  COEFFICIENT 

Aerosol  extinction  includes  absorption  and  scattering 
by  aerosol  particles.   Although  scattering  by  aerosol 
particles  whose  radius  is  smaller  than  about  0.03  times  the 
wavelength  of  the  light  may  be  calculated  by  Rayleigh 
theory,  for  the  whole  range  of  particle  radius  Mie  theory 
must  be  applied.   The  bases  for  Mie  scattering  are  the 
interactions  between  electromagnetic  waves  and  the  elec- 
tric charges  that  constitute  matter.   In  the  molecular 
scattering  case  which  was  described  before,  only  a  single 
dipole  is  involved,  but  in  the  case  of  aerosol  scattering, 
a  particle  consists  of  many  closely  packed,  complex  mole- 
cules.  Thus  it  may  be  considered  an  array  of  multipoles, 
since  an  aerosol  particle  is  a  dispersed  system  of  small 
particles  suspended  in  the  air,  such  as  combustion  pro- 
ducts, dust  grains,  bits  of  sea  salt  surrounded  by  water 
molecules,  volcanic  ash,  etc.   The  multipoles  give  rise  to 
secondary  electric  and  magnetic  partial  waves.   In  the  far 
field  two  kinds  of  waves  combine  to  produce  the  scattered 
wave.   Interference  between  two  waves  is  mainly  determined 
by  the  wavelength  of  incident  light,  and  the  size  and  rela- 
tive refractive  index  of  the  particle.   A  useful  parameter 
used  in  Mie  theory  is  the  size  parameter  a,  a  dimension- 
less  quantity  defined  by  the  relation 


a      =   2ttt   =   kr  (2-27) 
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where  r  is  the  radius  of  the  scatterer  or  absorber,  and  X 
is  the  wavelength  of  the  incident  radiation.  a    is  the 
parameter  on  which  aerosol  scattering  and  absorption  depend 
most  strongly,  rather  than  on  radius  of  particle  alone  or 
on  the  wavelength  of  incident  radiation  alone.   The  refrac- 
tive indices  used  in  calculations  are  the  relative  refrac- 
tive indices  of  the  aerosol  particles  with  respect  to  the 
refractive  index  of  the  medium  surrounding  the  particles. 
Since  air  has  refractive  index  very  close  to  unity,  the 
relative  refractive  index  is  assumed  equal  to  the  absolute 
refractive  index  of  the  aerosol  particle,  and  is  generally 
represented  by  a  complex  number,  where  the  imaginary  part 
takes  care  of  possible  absorption  by  the  particle.   The 
refractive  indices  of  different  kinds  of  particles  used  in 
calculations  are  reproduced  and  presented  in  table  10 
[Ref .  12] ,  taking  into  account  the  wavelength  dependence. 
Total  light  scattered  or  scattered  and  absorbed  by  a 
single  particle  is  given  by  its  total  scattering  or 
extinction  cross  section  which  is  related  to  the  geometric 
cross  section  by  a  parameter  generally  called  the  effi- 
ciency factor.   According  to  the  Mie  theory,  the  scattering 
and  extinction  efficiency  factors  are  given  by  the  equations 


q    -  !5=   =  \      l    (2n+l)Man!2+|bn|2]     (2-2S) 

01   n=l 


ex  "  ^7  =  -T  l    (2n+l)[Re(an+bn)l       (2-29 
"r       a   n=l 


Thus  the  efficiency  factor  is  defined  as  the  ratio  of 
apparent  cross  section  to  geometric  cross  section  of  the 
particle.   asc,  a^  in  the  above  are  the  apparent  scattering 
and  extinction  total  cross  sections  respectively.   Re 
represents  the  real  part  of  the  argument.   The  coefficients 
an#  b   are  the  complex  functions  which  represent  the 
amplitude  of  n    electric  partial  wave  and  n   magnetic 
partial  wave  respectively.   These  coefficients  may  be 
found  from  Ricatti-Bessel  functions  for  a  given  complex 
refractive  index  and  size  parameter  a. 

The  scattering  efficiency  factor  approaches  the  value 
of  two  in  the  short  wavelength  approximation,  which  implies 
that  the  scattering  cross  section  is  almost  twice  as  large 
as  the  geometric  cross  section  of  the  particle. 

Since  the  infinite  series  represented  by  equations 
(2-28)  and  (2-29)  are  slowly  converging,  calculations  are 
generally  computer  time  consuming.   For  particles  with  real 
part  of  refractive  index  between  1  and  1.5  and  imaginary 
part  of  refractive  index  between  0  and  0.25,  the  Van  de 
Hulst  approximation  may  be  used  with  some  correction  fac- 
tors.  Within  the  approximation  range,  it  was  reported  by 
Deirmendjian  [Ref.  13]  that  the  efficiency  factor  could 
be  predicted  within  the  precision  range  of  ±  0.05  Q.  • 

GX 

The  method  of  calculation  using  the  Van  de  Hulst 
approximation  is  presented  in  Appendix  C. 

The  efficiency  factor  of  a  single  particle  then  may 
be  extended  for  use  in  calculation  of  the  aerosol  extinction 
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coefficient  due  to  many  particles  for  a  given  size 
distribution  such  that 


a2 


3a   =   —   /a   Qsc  n(a)  da  (2-30) 

K 

a2 
V,      =  -\   /a2  Q   n(a)  da  (2-31) 

al 


where  K  is  the  propagation  constant  and  n(a)  is  the  size 

distribution  defined  by  equations  (1-6)  and  (2-27) . 

Once  the  aerosol  scattering  and  extinction  coefficients 

are  determined,  then  the  aerosol  absorption  coefficient  may 

be  found  from  K   =  y   -  3  ,  if  needed. 

a    a    a 

Size  distributions  used  here  for  calculation  of  aerosol 
extinction  coefficient  are  a)  equation  (l-6a)  which  is 
independent  of  wind  velocity  and  relative  humidity  having 
75%  sea  spray  and  25%  rural  aerosol  mixture  which  is  com- 
posed of  30%  dustlike  and  70%  water  soluble  substance, 
b)  equation  (l-6b)  which  adds  the  effects  of  wind  velocity 
and  relative  humidity. 

If  the  total  number  of  particles  is  not  known,  then  the 
size  distribution  n(a)  or  n(r)  is  generally  normalized  to 
one  particle  per  cubic  centimeter.    Whatever  value  the 
size  distribution  is  normalized  to  we  can  put  a  normaliza- 
tion constant  into  equations  (2-30)  and  (2-31) .   Thus  we 
may  write 
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a2 


S    =   G  -~   /   a   Q    n(a)  da 
a       vJ  sc 


al 


where  C  is  the  normalization  constant  and  may  be  found 
from  the  relation  [Ref.  19] 


R   =      3"9 


8.(0.55  urn) 
a 


Here  R  is  the  visibility  in  kilometers  providing  8  (0.55) 

a 

in  inverse  kilometers. 

Therefore  knowing  the  range  R  for  visibility,  the  aero- 
sol scattering  coefficient  can  be  found  from  the  equation 
above.   The  aerosol  scattering  coefficient  at  0.55  micro- 
meters calculated  by  equation  (2-30)  then  allows  determination 
of  the  value  of  the  normalization  constant  C.   Once  C  is 
determined  at  0.55  micrometers,  it  is  the  same  for  the 
whole  range  of  wavelengths. 

Since  aerosol  absorption  by  particles  is  much  less  than 
scattering,  the  visibility  range  may  also  be  written  as 


R   =      3"9 


U   (0.55  ym) 
a 


The  curves  of  extinction  coefficient  versus  wavelength 
for  varied  wind  speed  and  relative  humidity  shown  as  figures 
10  and  11  are  normalized  so  that  the  total  aerosol  extinction 
coefficient  at  0.55  ym  is  0.17  km   .   This  corresponds  to 
a  visibility  of  23  kilometers. 
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III.   ALGORITHM  FOR  CALCULATION  OF 

MOLECULAR  ABSORPTION  COEFFICIENT 

The  computer  program  presented  as  Appendix  A  is 

intended  to  calculate  the  molecular  absorption  coefficient 

including  the  water  vapor  continuum  effect.   As  described 

in  Section  II. B,  the  program  reads  the  line  parameters 

(v .  . ,    a .  . ,    S..)  from  a  magnetic  taoe.   However,  other 
*  l]    lj    l]  3 

inputs  of  the  program  are  the  amounts  of  the  air  molecules 
which  have  to  be  specified  in  the  units  of  molecules  per 
square  centimeter  per  kilometer  path,  and  they  should  be 
supplied  by  the  user. 

The  amounts  of  different  types  of  air  molecules  excluding 
the  ozone  and  water  vapor  are  calculated  from  the  ideal  gas 
law, 

P   =   p  R  T 

where  P  is  air  pressure  in  Newtons  per  square  meter,  T 
is  air  temperature  in  degrees  Kelvin,  R  is  the  gas  con- 
stant for  air  in  Joules  per  kilogram  per  degree  Kelvin 
and  has  the  value  of  287.06  Joules  per  kilogram  per  degree 
Kelvin.  For  given  values  of  pressure  and  temperature,  the 
density  of  air  can  be  found  in  the  units  of  kilograms  per 
cubic  meter.   Then  the  density  of  air  is  used  to  find  the 
amounts  of  the  uniformly  mixed  gases  (02,  N2 ,  C02 ,  CO, 
N20,  CHJ  such  that 
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P.   =   X.  10~6  p  . 
i       i       air 


where  Xi  is  the  concentration  in  PPM  by  volume  taken  from 
table  4,  and  p^  is  the  density  or  the  mass  per  unit  volume 
of  the  corresponding  constituent.   After  finding  the  amount 
of  a  given  type  of  molecule  in  air  in  the  units  of  mass 
per  unit  volume,  it  is  necessary  to  convert  it  to  the 
units  of  molecules  per  unit  volume  to  be  able  to  use  it 
in  the  computer  program.   This  can  be  done  through  the 
relation 


P  •  N,. 

l   A 

W.   = 


i  -3  ' 

M.  10 

l 


N. 


where  p.  has  the  units  of  kilograms  per  cubic  meter, 

+  23 

is  Avagadro ' s  number  which  is  6.02  3  x  10     molecules  per 

mole  and  M.  is  the  molecular  weight  in  grams  per  mole  of 
the  corresponding  constituent.   Now  the  amount  of  the  i 
constituent  of  air  is  given  by  the  formula  above  in  units 
of  molecules  per  cubic  meter.   Final  conversion  may  be 
made  into  units  of  molecules  per  square  centimeter  per 
kilometer  which  may  be  used  directly  as  input  to  the  com- 
puter program  presented  in  Appendix  A.   W.  may  also  be 
expressed  in  terms  of  P  and  T  as 


19  P  Xi 
W.   =   2.0982  x  iox:7  £  -±  3-1, 

1  T  M. 
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where  P  is  air  pressure  in  millibars,  T  is  temperature  in 

degrees  Kelvin,  X.  and  M.  are  as  defined  before.   Now  W. 

1       r  1 

is  the  amount  in  molecules  per  square  centimeter  per 
kilometer  for  the  i    type  of  air  molecules. 

Using  the  equation  (3-1) ,  the  amount  of  carbon  dioxide 
which  has  X±  =  330  and  M.  =  44  is  5.5943  *  1020  molecules 
per  square  centimeter  per  kilometer  in  the  air  where 
P  =  1017.5  millibars  and  T  =  286.22  degrees  Kelvin. 

The  densities  of  non  uniformly  mixed  gases  are  more 

complicated  to  determine.   For  ozone  at  sea  level,  it  is 

—  5  16 

taken  as  6  *  10    grams  per  cubic  meters  or  7.5288  x  10 

molecules  per  square  centimeter  per  kilometer  at  any  given 

temperature  and  pressure. 

The  amount  of  water  vapor  and  its  partial  pressure  can 
be  found  from  the  dew  point  temperature  or  relative  humidity 
in  the  air  [Ref.  15].   Different  dew  point  temperatures 
and  corresponding  density  and  partial  pressure  of  water 
vapor  are  tabulated  in  table  11. 

The  water  vapor  density  determined  by  interpolating 
from  table  11  may  be  converted  into  number  density  (per 
square  centimeter  per  kilometer  path)  as  previously 
described. 
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IV.   RESULTS 

Although  several  band  models  were  introduced  by  many 
workers  in  order  to  fit  the  molecular  absorption  lines  in 
the  atmosphere,  line  by  line  calculations  are  made  here, 
making  use  of  the  experimentally  found  and  compiled  data. 
The  molecular  absorption  coefficient  is  the  most  important 
attenuation  coefficient  to  be  calculated,  since  mainly  it 
is  the  one  which  determines  the  atmospheric  windows  which 
are  good  enough  for  propagation  of  light.   Also  it  depends 
strongly  on  the  frequency  of  propagation.   For  any  fre- 
quency of  laser  radiation,  the  molecular  absorption  coeffi- 
cient can  not  be  computed  by  interpolating  between  two  close, 
known  values.   It  can  be  computed  either  by  the  line  by 
line  method  (infinite  resolution)  or  by  interpolation  from 
calculations  with  resolution  less  than  0.01  cm   .   Therefore 
the  reported  values  of  molecular  absorption  coefficients 
may  be  used  directly.   If  any  other  frequency  or  condition 
is  needed,  the  computer  program  presented  as  Appendix  A  may 
be  used.   Atmospheric  windows  in  the  3-5  micrometer  and  8-14 
micrometer  regions  are  mainly  available  for  atmospheric 
transmission.   HF/DF  laser  wavelengths  fall  in  the  near 
and  middle  infrared  region  of  the  spectrum.   In  particular 
the  P2(8)  transition  of  the  DF  laser  with  wavelength  3.8007 
micrometers  shows  relatively  high  transmittance  falling 
in  the  3-5  micrometer  region  of  the  spectrum.   The  molecular 
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absorption  coefficient  for  this  particular  wavelength  is 
0.0104  km    in  the  midlatitude  winter  model  where  the 
molecular  scattering  coefficient  is  essentially  zero. 
This  shows  99  percent  transmittance  along  a  one  kilometer 
path. 

Calculations  intended  for  use  in  the  propagation  studies 
in  the  Monterey  Bay  area  make  use  of  the  known  range  of 
conditions  in  Monterey,  taking  into  account  the  yearly  and 
monthly  averages  of  air  pressure,  air  temperature  and  dew 
point  temperature  in  the  Monterey  Bay  area.   These  yearly 
and  monthly  averages  are  taken  from  [Ref.  16]  and  [Ref.  17]; 
the  quantities  used  in  calculations  are  tabulated  in  table 
12.   In  this  table  the  mean  temperature  is  used  which  is 
the  mean  value  of  maximum  and  minimum  average  temperatures 
corresponding  to  day/night  condition. 

Also  the  amounts  of  air  constituents  used  in  the  calcu- 
lations are  reported  in  table  17. 

Table  9  shows  the  calculated  molecular  absorption  for 
some  selected  wavelengths  for  the  midlatitude  winter  model, 
using  the  computer  program  presented  in  Appendix  A.   Table 
13  contains  the  molecular  absorption  coefficients  for  yearly 
averages  in  the  Monterey  Bay  area. 

Tables  14,  15,  16  show  the  molecular  scattering  and 
absorption  at  1.06,  3.8007  and  10.591033  micrometers 
respectively  for  the  monthly  averages  in  the  Bay  area. 
Since  the  molecular  scattering  is  negligible,  it  is  assumed 
to  be  zero  at  3.8007  and  10.591033  micrometers. 
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The  aerosol  extinction  coefficients  are  calculated 
for  the  rural  and  maritime  models  using  the  size  distri- 
bution given  by  equation  (l-6a)  and  presented  in  figure  3 
and  figure  9  respectively.   The  maritime  model  with  7  5% 
sea  spray  is  normalized  so  that  the  total  number  of  parti- 
cles is  about  44  per  cubic  meter.   These  results  agree  with 
the  previously  reported  values  [Ref.  12]  within  1.8  percent. 

The  relative  humidity  and  wind  velocity  dependent  size 
distribution  is  used  in  figure  10  and  figure  11,  normalizing 
the  curves  to  a  visibility  of  23  kilometers.   In  both 
figures  71%  sea  spray  and  29%  rural  aerosols  are  assumed. 
Although  at  0.55  micrometers,  the  rural  aerosol  refractive 
indices  are  outside  the  approximation  range,  the  results 
seem  to  be  valid  since  the  percentage  of  rural  aerosols 
is  low. 

As  seen  in  figure  10,  increasing  wind  speed  increases 
the  attenuation  in  the  longer  wavelength  region.   Similarly 
in  figure  11,  increasing  relative  humidity  increases  the 
attenuation  in  the  longer  wavelength  region. 

Figures  12  through  19  show  the  transmittance  along  a 
path  length  of  10  kilometers,  due  to  the  molecular  absorption 
only. 

Figures  12,  14,  16,  18  are  produced,  using  a  spectral 
resolution  of  0.01  cm    for  use  with  laser  sources. 

A  degraded  resolution  is  used  in  figures  13,  15,  17,  19 
by  averaging  the  transmittance  at  each  point  over  a  spectral 
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bandwidth  of  0.5  cm   .   Therefore  the  figures  13,  15,  17, 
19  are  for  use  with  the  broadband  sources  only. 

Each  one  of  the  figures  12,  13,  16,  17,  18,  19  covers 
a  spectral  interval  of  about  25  cm"1  near  the  wavelength  of 
3.8  micrometers  (i.e.,  2631. 09  cm"1) .   The  figures  14  and 
15  are  for  the  region  about  the  wavelength  of  10.6  micrometers 
(i.e. ,  943.936  cm"1) . 

Figures  12  through  15  show  the  transmittances  for 
yearly  averages  in  the  Monterey  Bay  area. 

Figures  16  through  19  are  produced  for  the  different 
partial  pressures  of  water  vapor,  having  the  same  air 
temperature  of  290  degrees  Kelvin  and  the  air  pressure  of 
1018  millibars.   The  partial  pressure  of  water  vapor  in 
figures  16  and  17  is  7.48  mbars,  and  it  is  15.48  mbars  in 
figures  18  and  19. 

Using  the  high  resolution  of  0.01  cm    and  a  degraded 
resolution,  a  sensitivity  analysis  due  to  the  air  pressure, 
temperature  and  partial  pressure  of  water  vapor  has  been 
carried  out  near  the  wavelength  of  3.8  micrometers. 

It  has  been  observed  that  the  change  in  the  air 
pressure  and/or  temperature  which  could  affect  the  line 
widths,  line  intensities  and  amounts  of  the  absorbing  air 
molecules,  does  not  contribute  to  the  transmittance  signi- 
ficantly, at  the  region  near  3.8  micrometers.   The  reason 
could  be  the  small  number  of  absorbing  lines  which  are 
mostly  due  to  water  vapor.   Comparing  figures  12,  16,  18  or 
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figures  13,  17,  19  it  can  easily  be  seen  that  the  change 
of  the  partial  pressure  of  water  vapor  effects  the  trans- 
mittance  significantly.   8  millibars  change  in  the  partial 
pressure  of  water  vapor  changes  the  transmittance  of  the 
order  of  6-7  percent  at  3.8  micrometers  region. 
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APPENDIX  C 

CALCULATION  OF  EXTINCTION  EFFICIENCY  BY  USING 
THE  VAN  DE  HULST  APPROXIMATION 


This  approximation  was  first  developed  by  Van  de  Hulst 
[Ref.  13]  following  the  Huygens  principle.  It  is  good  for 
spherical  particles  of  arbitrary  size  but  |m|  ■*•  1. 

The  refractive  index  of  the  particle  is  denoted  by  m 
which  has  real  and  imaginary  parts. 


m  =   n   -in. 
r     1 


Defining  the  normalized  size  parameter  p  and  absorption 
parameter  g  such  that 


p   =   2a  (n   -  I! 


,   -i    ni     . 


the  extinction  efficiency  may  be  written  as 


Q'    =   2  -  4  cos  g  exp[-ptan  g]  sin  (p-g) 
ex  p 


+  4  (cos  9) 2{cos  2g  -  exp[-ptan  g]  cos(p-2g)} 


This  uncorrected  extinction  efficiency  Q'   overestimates 
the  cross  section  for  small  p  and  underestimates  it  in 
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varying  degrees  as  p  approaches  and  surpasses  the  value 

4.08,  corresponding  to  the  first  maximum  in  Q1  . 

ex 

A  correction  to  the  extinction  efficiency  is  made  in 
the  form 


Q     =   (1  +  D)  Q' 

ex  '    vex 


where  D,  the  correction  factor,  is  given  by  the  expressions 

(n  -1) 2  0.2p-n  +1 

D   =   1.632  nr[f(g)+1]  +  (nr-l)f(g)    when   P  <  5 (nr-l) 


D   "   8^T5-n-[f(9)+1^     ^en   5(nr-l)  <  p  <   j£g*   g 

r  a 


(nr  1}        .^..  .       4.08    „    ,      4.08 

D   "  -,„    TTTu^    -,   [f  (g)+U    when   r^rrr  -  <  P  < 


2n  (l+3tan  g)  L  v^'   J  l+3tan  g  -  '  -  1+tan  g 


2.04(n  -1)  Q 

D   =   ^4 [f(g)+l]     when  p  >  ,11     

n  pf (g)       ^   -  K    1+tan  g 


where 

2 


f(g)   =   1+4  tan  g  +  3 (tan  g) 


Within  the  limits  of 


1  <  n   <_  1.5 


0  <  n.  <  0. 25 

—   i  — 


66 


Q    is  accurate  within  an  error  of  about  ±  0.05  Q 
vex  ex 

for  a  wide  range  of  sizes  and  types  of  index  of  refraction 
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TABLE  1.   MODIFICATION  FACTOR  FOR  C02  LINES 


v  -  v  x 

o 

1/CM 


0.0  1,00 

0.5  1.00 

0.6  0.96 

0.7  0.89 

0.8  0.32 

0.9  0.77 

1.0  0.70 

1.2  0.60 

1.5  0.50 

2.0  0.41 

2.5  0.34 

3.0  0.31 

5.0  0.29 

3.0  0.23 

10.0  0.19 

15.0  0.00 
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TABLE    2 .       VIBRATIONAL    PARTITION    FUNCTIONS 


MOLECULE       J  TEMPERATURE    (DEGREES    KELVIN) 

175  200  225  250  275  296  325 


H20  1.5    1.0000  1.0000  1.0000  1.0000    1.0000    1.000  1.0000 

C02  1.0    1.0095  1.0192  1.0327  1.0502    1.0719    1.093  1.1269 

03  1.5    1.0040  1.0070  1.0130  1.0220    1.0330    1.046  1.0660 

N20  1.0    1.0170  1.0300  1.0480  1.0720    1.1000    1.127  1.1700 

CO  1.0    1.0000  1.0000  1.0000  1.0000    1.0000    1.000  1.0000 

CH4  1.5    1.0000  1.0000  1.0010  1.0020    1.0040    1.0C7  1.0110 

02  l.C   1.0000  1.0000  1.0000  1.0000    1.0000    1.000  1.0010 
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TABLE  3.   SELF  BROADENING  ABSORPTION  COEFFICIENTS  FOR 
WATER  VAPOR  IN  THE  3.5-4.2  MICROMETER  REGION 


Cs(v,296) 


1/CM         l/(PR.CM)ATM. 


2350  0.230 

2400  0,187 

2450  0.147 

2500  0.117 

2550  0.097 

2600  0.087 

2650  0.100 

2700  0.120 

2750  0.147 

2800  0.174 

2850  0.200 

2900  0.240 

2950  0.230 

3000  0.330 
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TABLE  4.   CONCENTRATIONS  OF  UNIFORMLY  MIXED  GASES 
IN  THE  AIR 


CONSTITUENT     CONCENTRATION 

P°M  BY  VOLUME 


C32  3.300E+2 

N23  2.300E-1 

CD  7.500E-2 

CH4  1.600E+0 

02  2.095E+5 

N2  7.808E+5 
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TABLE  5.   MODELS  OF  ATMOSPHERE  AT  SEA  LEVEL 


MODEL  PRES.  TEMP.  DENSITY      H20  03 

MB.  DEG.    KEL.     GR/M3         GR/M3       GR/M3 


TROPICAL  1013  300.0  1167  19.0  5.6E-5 


MIDLATITUDE  1013  294.0  1191  14.0  6.0E-5 

SUMMER 

MIDLATITUDE  1013  272.2  1301  3.5         6.0E-5 

WINTER 

SUBARCTIC  1010         287.0  1220  9.1  4.9E-5 

SUMMER 

SUBARCTIC  1013         257.1  1372  1.2         4.1E-5 

WINTER 

U.S.    STANDARD  1013  283.1  1225  5.9         5.4E-5 

ATM.    ,1962 
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TABLE  6.   SIZE  DISTRIBUTION  PARAMETERS  OF  AEROSOL  MODELS 


MODEL  Nx  rx  ax  N2  ^2 


RURAL  0.9999975  0.005u  0.475       2.5E-6       0.5'M        0.475 

URBAN  0.9999975  0.005U  0.475       2.5E-6       0  .  5P        0.475 

MARITIME  1  0.  005U  0.475 

(CONTINENTAL) 

MARITIME  1  0.3y  0.4 

(SEA    SPRAY) 


N       IS    NORMALIZED    TO    1    PARTICLE/CM3 
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TABLE  7.   HF  LASER  FREQUENCIES   [10] 

VIBRATIONAL           ROTATIONAL  EMISSION 
TRANSITION    ID         TRANSIT ICN    10       FREQUENCY    (1/CM) 

1-0                                  P6  3693.50 

P7  3644.16 

P8  3593.80 

P9  3542.20 

P10  3489.59 

Pll  3436.12 

P12  3381.50 

P13  3326.21 

P14  3269.90 

P15  3212.80 

2-1                                P2  3708.86 

P3  3666.35 

P4  3622.71 

P5  3577.47 

P6  3531.31 

P7  3483.63 

P8  2425.17 

P9  3335.34 

P10  3334.55 

Pll  3282.36 

P12  3230.13 

P13  3176.60 

P14  3122.14 

P15  3067.22 

3-2  '                             P2  3544.51 

P3  3503.80 

P4  3461.54 

P5  3413.16 

P6  3373.46 

P7  3327.73 

P3  3230.64 

4-3                               P5  3262.53 

P6  3219.50 
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TABLE    7.       HF    LASER   FREQUENCIES     [10]     (Continued) 

P7  3175.34 

P8  3130.09 

P9  3083.33 

5-4  P4  3150.67 

P5  3110.34 

P6  3068.63 

P7  3026.21 

P8  2982.51 

P9  2937,79 

6-5  P5  2961.68 

P6  2921.74 

P7  2380.70 

P3  2838.59 
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TABLE    8.       DF    LASER   FREQUENCIES     [4] 

VIERATIONAL  ROTATIONAL  EMISSION 

TRANSITION       ID      TRANSITICN    ID      FREQUENCY    (1/CM) 

1-0  PI  2384. 934 

P2  2362.652 

P3  2839.779 

P4  2816.362 

?5  2792.437 

P6  2767.914 

P7  2743.028 

P8  2717.536 

P9  2691.409 

P10  2o65.200 

Pll  2633.396 

Pll  2611.125 

P10  2534.910 

P14  2557.090 

P15  2527.060 

P16  2498.020 

2-1  P3  2750.050 

P4  2727.330 

P5  27C3.980 

P6  2680.280 

P7  2655.970 

P8  2631.090 

P9  2605.370 

P10  2530.160 

Pll  2553.970 

Pll  2527.470 

P13  2500.320 

P16  2417.270 

2-2  P3  2662.170 

P4  2640.040 

P5  2617.410 

P6  2594.230 

P7  2570.510 
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TABLE    8.       DF   LASER   FREQUENCIES     [4]     (Continued) 

P8  2546.370 

P9  2521.810 

P10  2496.610 

Pll  2471.340 

P12  2445.290 

P13  2419.020 

P14  2392.460 

4-3  P5  2532.500 

P6  2509.860 

P7  2486.830 

P8  2463.250 

P9  2439*290 

PIC  2414.890 

5-4  P7  2404.630 

7-6  P8  2222.680 

P10  2177.990 

Pll  2155.030 

P12  2131.680 

8-7  P7  2165.930 

P8  2144.300 

P9  2123.240 

P10  2101.270 

Pll  2056.140 

P10  2033.010 

9-8  P6  2108.480 

P7  2038.340 

P8  2067.760 

P10  2025.360 

Pll  2003.560 

P12  1981.330 
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TABLE  9.   MIDLATITUDE  WINTER  MODEL  MOLECULAR  ABSORPTION 
COEFFICIENTS  FOR  INFINITE  RESOLUTION 


FREQUENCY  MOLECULAR    ABSORPTION 

(1/CM)  COEFFICIENT     Q/KM) 

924.975  4.27E-1 

942.384  6.17E-2 

944.195  6.6IE-2 

945.981  6.665-2 

9  47.743  7.22E-2 

959.393  3.45E-2 

973.289  7.49E-2 

974.623  6.96E-2 

975.931  1.32E-1 

977.215  7.28E-2 

1046.854  8.27E-2 

1048.661  8.97E-2 

1050.441  9.81E-2 

1075.988  3.90E-2 

1077.303  8.21E-2 

1078.591  7.37E-2 

2177.990  3.39E-2 

2414.890  1.16E-2 

2463.250  1.26E-2 

2611.125  7.19E-2 

2617.410  7.19E-3 

2631.090  1.04E-2 

2638.396  7.90E-2 

2640.040  1.67E-2 

2655.970  2.24E-2 

2662.170  1.24E-2 

2665.200  1  .64E-2 

2703.980  9.55E-3 

2717.536  3.37E-2 

2727.380  1.38E-2 

2743.023  1.16E-2 

2750.050  1.61E-2 


78 


TABLE    10.       REFRACTIVE    INDICES    OF   AEROSOL    PARTICLES     [Ref.     12] 

WAVELENGTH    WATER    SOLUBLE    DUST    LIKE         SQCT    LIKE  OCEANIC 

(pm)  nr  ni  nr  ni  nr        ni        n  n. 


0.2  1.530  0.070  1.53    0.070  1.50    0.35  1.429    2.87E-5 

0.25  1.530  0.030  1.53    0.030  1.62    C.45  1.404    1.45E-6 

0.3  1.530  0.008  1.53    0.008  1.74    0.47  1.395    5.83E-7 

0.337.:,  1.530  0.005  1.53   0.008  1.75    C.47  1.392    1.20E-7 

0.4  1.530  0.005  1.53    0.003  1,75    C.46  1.385    9.90E-9 

0.48S  1.530  0.005  1.53    0.003  1.75    0.^5  1.382    6.41E-9 

0.5145  1.530  0.005  1.53    0.003  1.75    0.45  1.381    3.70E-9 

0.55  1.530  0.006  1.53    0.003  1.75    0.44  1.331    4.26E-9 

0.6323  1.530  0.006  1.53    0.003  1.75    0.43  1.377   1.52E-3 

0.6943  1.530  0.007  1.33    0.008  1.75    0.43  1.376    5.04E-3 

0.36  1.520  0.012  1.52    0.008  1.75    0.43  1.372    1.09E-6 

1.06  1.520  0.017  1.52  0.003  1.75  C.44  1.367  6.01E-5 
1.3  1.510  0.020  1.46  0.008  1.76  C.45  1.365  1.41E-4 
1.536  1.510  0.023  1.40  0.008  1.77  0.46  1.359  2.43E-4 
1.3  1.460  0.017  1.33  0.008  1.79  0.43  1.351  3.11E-4 
2.0  1.420  0.003  1.26  0.003  1.80  0.49  1.347  1.C7E-3 
2.25  1.420  0.010  1.22  0.009  1.31  0.50  1.324  5.50E-4 
2.5  1.420  0.012  1.13    0.009  1.82    0.51  1.309    2.39E-3 

2.7  1.400  0.055  1.13  0.013  1.83  0.52  1.249  1.56E-2 
3.0  1.420  0.022  1.16  0.012  1.34  C.54  1.429  1.97E-1 
3.2  1.430  0.008  1.22  0.010  1.36  C.54  1.481  6.69E-2 
3.3923  1.430  0.007  1.26  0.013  1.37  0.55  1.439  1.51E-2 
3.5  1.450  0.005  1.2  3  0.011  1.38  0.56  1.423  7.17E-3 
3.75  1.452  0.004  1.27  0.011  1.90  0.57  1.398  2.90E-3 
4.0  1.455  0.005  1.26  0.012  1.92  0.58  1.388  3.69E-3 
4-.S  1.460  0.013  1.26  0.014  1.94  0.59  1.377  9.97E-3 
5.0  1.450  0.012  1.25  0.016  1.97  0.60  1.366  9.57E-3 
5.5  1.440  0.013  1.22  0.021  -.99  0.61  1.333  9.31E-3 
6.0  1.410  0.022  1.15  0.037  2.02  0.62  1.306  7.96E-2 
6.2  1.430  0.027  1.14  0.039  2.03  0.63  1.421  6.91E-2 
6.5  1.460  0.033  1.13  0.042  2.04  0.63  1.374  2.94E-2 
7.2  1.400  0.070  1.4C    0.055  2.06    0.95  1.543    2.49E-2 
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TABLE    10.       REFRACTIVE    INDICES    OF   AEROSOL    PARTICLES     (Continued) 
[Ref.     12] 

7.9  I. 200    C.065  1.15  0.040    2.12  0.67  1.324    2.79E-2 

3.2  1.010    0.100  1.13  0.074    2.13  0.63  1.324    3.086-2 

3.5  1.300    0.215  1.30  0.090    2.15  C.69  1.336    3.36E-2 

8.7  2.400  0.290  1.40  0.100  2.16  C.69  1.336  3.56E-2 
9.0  2,560  0.370  1.70  0.140  2.17  0.70  1.373  3.65E-2 
9.2  2.200  0.420  1.72  0.150  2.13  0.70  1.356  3.71E-2 
9.5  1.950    0.160  1.73  0.162    2.19  0.71  1.339    3.63E-2 

9.8  1.370  0.095  1.74  0.162  2.20  0.72  1.324  3.33E-2 
10.0  1.820  0.090  1.75  0.162  2.21  0.?2  1.310  4.06E-2 
10.591  1.760  0.070  1.62  0.120  2.22  0.73  1.271  5.22E-2 
11. 0  1.^20  0.050  1.62  0.105  2.23  C.73  1.246  7.31E-2 
11.5  1.670  0.047  1.59  0.100  2.24  C . 74  1.227  1.05E-1 
12.5  1.620  0.053  1.51  0.090  2.27  0.75  1.208  1.90E-1 
13.0  1.620  0.055  1.47  0.100  2.23  0.76  1.221  2.23E-1 
14.0  1.560    0.073  1.52  0.085    2.31  0.73  1.267    2.71E-1 

14.8  1.440  0.100  1.57  0.100  2.33  0.79  .1.307  2.92E-1 
15.0  1.420    0.200  1.57  0.100    2.33  0.79  1.321    2.97E-1 

16.4  1.750    0.160  1.60  0.100    2.36  0.81  1.407    3.31E-1 

17.2  2.030  0.240  1.63  0.100  2.33  C.32  1.487  3.41E-1 
18.0  1.980    0.180  1.64  0.115    2.40  G.83  1.525    3.41E-1 

13.5  1.850  0.170  1.64  0.120  2.41  0.33  1.526  3.39E-1 
20.0  2.120    0.220  1.63  0.220    2.45  0.35  1.560    3.24E-1 

21.3  2.060  0.230  1.77  0.230  2.46  C.36  1.563  3.13E-1 
22.5  2.000  0.240  1.90  0.230  2.48  0.37  1.579  3.16E-1 
25.0  1.880    G.280  1.97  0.240    2.51  0.39  1.596    3.13E-1 

27.9  1.840  0.290  1.39  0.320  2.54  0.91  1.612  3.20E-1 
30.0  1.320  0.300  1.30  0.420  2.57  C.93  1.614  3. 20E- 1 
25.0  1.920  0.400  1.90  0.500  2.63  C.97  1.597  3.16E-1 
40.0  1.360    0.50  0  2.10  0.600   2.69  1.00  1.582    5.61E-1 
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TABLE    11.       WATER   VAPOR    PRESSURE    AND    DENSITY    FOR    DIFFERENT 
DEW    POINT    TEMPERATURES 


DEW     POINT    TEMP.  VAPOR    PRESSURE  VAPOR    DENSITY 

(DEC    KEL.)  (MILIBARS)  (GRAMS/M3) 


233  0.189  0,176 

238  0.314  0.286 

243  0.509  0.453 

248  0.307  0.705 

253  1.254  1.074 

258  1.912  1.605 

263  2.863  2.358 

268  4.215  3.407 

273  6.108  4.347 

278  8.719  6.797 

233  12.273  9.399 

288  17.044  12.330 

293  23.373  17.500 

298  31.671  23.050 

303  42.430  30.380 

308  56.236  39.630 

315  73.777  51.190 
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TABLE  13.   MOLECULAR  ABSORPTION  AND  SCATTERING  COEFFICIENTS 
FOR  YEARLY  AVERAGES  IN  THE  MONTEREY  BAY  AREA 


FREQUENCY 

TRANSITION 

MOL.    ABS . 

MOL.    SCAT. 

TRANS.    A 

ID 

COEF. 

COEF. 

1    KM.     P 

tl/CM) 

(1/KM) 

(1/KM) 

m 

944.195 

P(20) 

CO  2 

0.1520 

— 

35.36 

942.3  84 

P(22) 

CO  2 

0.1460 

- 

86.44 

2471.340 

P(ll) 

OF 

0.0230 

- 

97.^2 

2496.610 

P(10) 

DF 

0.01  34 

- 

93.17 

2527.470 

P(12) 

OF 

0.  0172 

- 

93.29 

2553.970 

PUi) 

Dc 

0.  0199 

- 

93.03 

2530.160 

P  ( 1  0 ) 

DF 

0.02  93 

- 

97.11 

2605.370 

P(9  ) 

DF 

0.0340 

- 

96.65 

2631.090 

P(  3) 

DF 

0.0218 

- 

97.34 

2638.396 

PUI) 

DF 

0.2210 

- 

80.17 

2655. 97C 

P<7) 

DF 

0.0531 

- 

94.33 

2665. 20C 

PUO) 

DF 

0.0357 

- 

96.49 

26SC.28C 

P(6  ) 

DF 

0. 0547 

- 

94.63 

2703.980 

P(5) 

DF 

0.0208 

- 

97.94 

2717.536 

P(3) 

DF 

0.0930 

- 

91.12 

2743.023 

P(7  ) 

DF 

0.0  246 

- 

97.57 

2767.914 

P(  6  ) 

DF 

0.0594 

- 

94.23 

2792.437 

P<  5) 

OP 

0.0439 

- 

95.70 

3373.460 

P(6) 

HF 

0.1240 

- 

88.37 

343  5.1  7  C 

P(3  ) 

HF 

0.3470 

- 

70.71 

3483. 63 C 

P(7) 

HF 

2.1300 

- 

11.32 

3531.310 

P(6) 

HF 

86.5000 

- 

0.00 

3577.470 

P(5  ) 

HF 

35.0000 

- 

0.00 

3593  .300 

P(3  ) 

HF 

67.2000 

- 

0.00 

3622.710 

P(4) 

HJ= 

35.6000 

- 

0.00 

9433.96  2 

0.000014 

0.0003283 

99.92 

OATH 
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TABLE  14.   MOLECULAR  ABSORPTION  AND  SCATTERING  COEFFICIENT 
AT  1.06  MICROMETERS  IN  THE  MONTEREY  BAY  AREA 


MOL.     ABS.  MOL.     SCAT.       TOTAL    MOL.     TRANS.    4LCNG 

COEF.  COEF.  EXTINCTION!        1    KM.     PATH 

(1/KM)  (1/KMJ  (1/KM)  (%1 


JANUARY 

9  .270E-6 

9.39SE-4 

8.  486E-4 

99.92 

FEBRUARY 

9.700E-6 

8.377E-4 

8.474E-4 

99.92 

MARCH 

9.710E-6 

8.  3  49E-4 

3.4465-4 

99.<=2 

APRIL 

1 0.900E-6 

3.328E-4 

8.437E-4 

99.92 

MAY 

11 .800E-6 

3.272E-4 

S.390E-4 

99.92 

JUNE 

13.300E-6 

8.223E-4 

3.361E-4 

99.92 

JULY 

13.900E-5 

3.232E-i 

8.371E-4 

99.92 

AUGUST 

14.500E-6 

9.  220  E- 4 

8.365E-4 

99.92 

SEPTEMBER 

14.100E-6 

8.172E-4 

8.313E-4 

99.92 

OCTOBER 

12.500E-6 

8.236E-4 

3.3615-4 

99.92 

NOVEMBER 

10.400 E-6 

8.3  0^E-4 

8.403E-4 

99.92 

DECEMBER 

9  .410E-6 

8.344E-4 

3.  433E-4 

99.92 

YEARLY 

11.4005-6 

3.283E-4 

3.397E-4 

99.92 
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TABLE    15.       MOLECULAR   ABSORPTION    COEFFICIENTS    AT    3.8007 
MICROMETERS    IN    THE    MONTEREY    BAY    AREA 


MOL.     ABS.       TRANS.     ALONG      ORDER    OF 
COEF.  1    KM.    PATH       BEST    MONTH 

(1/KM)  {%)               TO    PROPAGATE 

JANUARY                      0.0194  93.08  2 

FEBRUARY                    0.0199  93.02  4 

MARCH                            0.0199  93.03  3 

APRIL                            0.0215  97.38  6 

MAY                                 0.0223  97.79  7 

JUNE                              0.0242  97.61  9 

JULY                              0.0249  97.54  11 

AUGUST                         0.0256  97.43  12 

SEPTEMBER                 0.0247  97.56  10 

0CT03ER                      0.0227  97.75  8 

NOVEMBER                    0.0202  93.00  5 

DECEMBER                    0.0193  98.09  1 

YEARLY                         0.0218  93.34 
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TABLE  16.   MOLECULAR  ABSORPTION  COEFFICIENTS  AT  10.591033 
MICROMETERS  IN  THE  MONTEREY  BAY  AREA 


TRANS.     ALONG      ORDER    OF 

1    KM.     PATH  BEST    MONTH 
(%)                TO    PROPAGATE 

JANUARY                      0.122                    33.51  2 

FEBRUARY                    0.123                    33.  01  4 

MARCH                            0.128                   33.02  3 

APRIL                           0.14-7                    36.37  6 


MOL. 

,    ABS. 

COEF. 

U/KM) 

0, 

.122 

0, 

123 

0, 

,12  8 

0. 

,14-7 

0. 

159 

0. 

,135 

0, 

,195 

0. 

204- 

0, 

,193 

0, 

,165 

0, 

133 

0, 

,121 

0. 

,152 

MAY  0.159  85.29  7 

JUNE  0.135  83.15  9 

JULY  0.19  5  32.32  11 

AUGUST  0.204  31.54  12 

SEPTEMBER  0.193  32.43  10 

OCTOBER  0.165  34.78  8 

NOVEMBER  0.133  37.53  5 

DECEMBER  0.121  38.57  1 

YEARLY  0.152  35.36 
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TABLE    18.       AEROSOL    EXTINCTION    COEFFICIENTS    FOR    MARITIME 

MODEL    WITH    7  5%    SEA    SPRAY    AND    25%    RURAL    AEROSOLS 


WAVELENJ7H  AEROSOL    EXTINCTION  AEROSOL    EXTINCTION 

(MICROMETERS)  COEFFICIENT    (1/KMJ  COEFFICIENT    (1/KM) 

CALCULATED  REPORTED    [Ref.     12] 


2.0  0.11330  0.11530 

2.5  0.09895  0.09962 

2.7  0.08253 

3.0  0.10430  0.10426 

3.2  0.10800 

3.3923  0.10250 

3.5  0.0^945  0.09899 

3.75  0.09299  0.09191 

4.C  0.03802  0.08670 

4.5  0.07875 

5.5  0.06007  0.05928 

6.0  0.05555  0.05485 

6.5  0.05825 

7.2  0.04304  0.04758 

7.9  0.04105  0.04063 

3.2  0.03958  0.03960 

8.5  0.03972  0.04045 
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FIGURE  1.   a.   MODEL  OF  AN  ELEMENTAL  SCATTERER 

b.   CREATION  OF  AN  INDUCED  DIPOLE  MOMENT 
BY  AN  ELECTRIC  FIELD 
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INCIDENT 


FIGURE  2.   THE  GEOMETRY  OF  SCATTERING  WITH  THE 
SCATTERER  AT  THE  ORIGIN 
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1  Vacuum  pump 

2  Charcoal  filter 

3  Scrubber 

4  Scrubber-vacuum  valve 

5  By-pass  valve 

6  Cavity  pressure  gauge 

7  Cavity  vacuum  valve 

8  Optical  resonator  mirrors 

9  Brewster  windows 

10  Cavity 

11  Electric  discharge  tube 

12  H2  or  D2  gas  tank 

13  SFg ,  02/  He  gas  tanks 

14  Ballast  resistors 

15  Ballast  resistor  cooling  fans 

16  D.C.  power  supply 

17  Water  tank 

18  Cooling  water  pump 

19  Oil  filter 

20  Gas  outlet  to  atmosphere 

FIGURE  5.  HF/DF  LASER  SYSTEM  BLOCK  DIAGRAM  (Continued) 
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